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Abstract. - We unify two paradigmatic mesoscopic mechanisms for the emergence of nonex- 
tensive statistics, namely the multiplicative noise mechanism leading to a linear Fokker-Planck 
(FP) equation with mhomogenous diffusion coefficient, and the non-Markovian process lead- 
ing to the nonlinear FP equation with homogeneous diffusion coefficient. More precisely, we 
consider the equation dp ^g t '^ = — -^\F(x)p{x,t)] + ^D-^ I [(j)(x,p)p(x,t)}, where D £ 1Z and 
F(x) = —dV(x)/dx, V(x) being the potential under which diffusion occurs. Our aim is to 
find whether <j>(x,p) exists such that the inhomogeneous linear and the homogeneous nonlinear 
FP equations become unified in such a way that the (ubiquitously observed) q-exponentials re- 
main as stationary solutions, ft turns out that such solutions indeed exist for a wide class of 
systems, namely when <fi(x,p) = [A + BV(x)] e [p(x, t)] v , where A, B, 6 and n are (real) constants. 
Our main result can be sumarized as follows: For 9/1 and arbitrary confining potential V(x), 



p(x,oo) « {l-0(l-q)V(x)} 



1/(1-9) 



-PV(x) 



where q = 1 + n/(9 — 1). The present approach 



unifies into a single mechanism, essentially long memory, results currently discussed and applied 
in the literature. 



One of the cornerstones of statistical mechanics is the 
functional connection of the thermodynamic entropy with 
the set of probabilities {pi} of microscopic configurations. 
[ For the celebrated Boltzmann-Gibbs (BG) theory, this 

'central functional is given by Sbg = — kJ2i=iPi^ n Pii 
where W is the total number of microscopic states which 
are compatible with the information that we have about 
the system. This powerful connection is in principle appli- 
cable to a vast class of relevant systems, including (clas- 
sical) dynamical ones whose maximal Lyapunov expo- 
nent is positive, thus generically warranting strong chaos, 
hence mixing in phase space, hence ergodicity (in some 
sense, Boltzmann embodied all these features in his in- 
sightful molecular chaos hypothesis). Within this the- 
ory, it ubiquitously emerges the Gaussian distribution 
Pg oc e~ l3x ((3 > 0). Indeed, this important proba- 
bilistic form (i) maximizes the (continuous version of the) 



entropy Sbg = —k J dxp (x) In [p (x)] under the basic con- 
straints of normalizability and finite width; (ii) constitutes 
the exact solution, for all values of space and time, of the 
simplest form of the (linear and homogeneous) Fokker- 
Planck equation, in turn based on the simplest form of the 
Langevin equation (which includes additive noise); (iii) 
is the N — > oo attractor of the (appropriately centered 
and scaled) sum of N independent (or weakly correlated 
in an appropriate sense) discrete or continuous random 
variables whose second moment is finite (Central Limit 
Theorem, CLT); (iv) is the velocity distribution (Maxwell 
distribution) of any classical many-body Hamiltonian sys- 
tem whose canonical (thermal equilibrium with a thermo- 
stat) partition function is finite, i.e., if the interactions be- 
tween its elements are sufficiently short-ranged, or inexis- 
tent. The simplest probabilistic model which realizes these 
paradigmatic features is a set of N independent equal bi- 
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nary random variables (each of them taking say the values 
and 1 with probability 1/2). The probability of having, 



for fixed N, n l's is given by 



>-JV 



Its limiting 



! (N-n)\ 

distribution is, after centering and scaling, a Gaussian (as 
first proved by de Moivre and Laplace), and its (extensive) 
entropy is the BG one, since Sbg(N) = 7Vfcln2. 

What happens with the above properties when the cor- 
relations between say the elements of a probabilistic model 
are strong enough (in the sense that they spread over all el- 
ements of the system)? There is in principle no reason for 
expecting the relevant limiting distribution to be a Gaus- 
sian, and the entropy which is extensive (i.e., S(N) oc N 
for N 3> 1) to be Sbg- The purpose of the present paper 
is to focus on such and related questions for a class of sys- 
tems which are ubiquitous in natural, artificial and even 
social systems, namely those which are scale-invariant in 
a probabilistic sense which we shall define below. Let us 
now discuss the frequent emergence of (/-exponentials, de- 
fined as 

P q (x)=N g [l-(l-q)pV(x)] 1 /^ =N q e-e v M , (1) 

where N q is a normalization factor; Pi(x) = N\e~P v ( x > 
is the standard BG case; for V(x) oc x 2 , P g (x) is a q- 
Gaussian, which displays asymptotic power-laws and can 
be seen as a natural generalization of the Gaussian (g> = 1). 

At this point let us make a few remarks. (i) q- 
Gaussians appear as the exact solutions of paradigmatic 
non-Markovian Langevin processes and their associated 
Fokker-Planck equations . Langevin equations with both 
additive and multiplicative noise [T], or Langevin equa- 
tions with long-range- memory [2], lead respectively to in- 
homogeneous linear [3] , or homogeneous nonlinear [4 , 5 
Fokker-Planck equations (see also [SHlOj). (ii) g-CLT at- 
tractors are g-Gaussians [llj . (iii) The extremization of 
the entropy S q with norm and finite width constraints 
yields q-Gaussians, where S q is a generalization of BG 
entropy, namely [T2"l[T3"] 



Sn — k 



1 - / dx\p(x)] q 
q~ i 



(q eR;S 1 = S BG ) (2) 



This entropy is, for q ^ 1, nonadditive (see [14] for the 
current definition of additivity), i.e., for arbitrary prob- 
abilistically independent systems A and B, the equality 
S(A + B) = S(A) + S(B) is not satisfied. However, for 
many systems a value of q, denoted by q en t, exists for 
which S qent is extensive, i.e., S qent (N) oc N (N ^> 1). 
As is well known, for all standard short-range-interacting 
many-body Hamiltonian systems, we have q cn t = 1. How- 
ever, some systems exist for which q cnt < 1 [L5HT8] . (iv) 
Numerical indications [19. for the distributions of veloci- 
ties in quasistationary states of long-range Hamiltonians 
[20] suggest g-Gaussians. Further, experimental and ob- 
servational evidence for g-Gaussians exists for the motion 
of biological cells [21H22] , defect turbulence [23] , solar wind 
[24 25 , cold atoms in dissipative optical lattices [2H] j dusty 



plasma [37] , among others (see also 28 ) . Numerical indi- 
cations are also available at the edge of chaos of unimodal 
maps [29] . 

A domain where the nonadditive entropy S q can be 
naturally incorporated is for describing anomalous diffu- 
sion like-phenomena. From modified Langevin equations, 
inhomogeneous linear or homogeneous nonlinear Fokker- 
Planck equations have been derived and used in order to 
obtain the mcsoscopic dynamic evolution of systems where 
such diffusion occurs. 

Let us consider here the nonlinear Fokker-Planck (FP) 
equation given by 



dp(x, t) d r „, . . 1 „ d 2 , , , 

' = - — [F(x)p(x,t)} + -D—tt(x,p)p(x,t)} 



dt 



dx 



x (3) 

where D e TZ is the coeficient of diffusion, F — — dv ^ the 

' ox 

drift term, and V a confining potential. We shall further 
assume the following wide connection: 



(x,p) = [g(V)} e \p(x,t)}\ 



(4) 



with g(V) = [A+BV], A, B, 9 and n being real constants. 

It should be noted that the present equation encom- 
passes also the general inhomogeneous (9 ^ 0) nonlin- 
ear (77 7^ 0) case. Our aim is to find whether a q- 
exponential P q {x) exists as a stationary solution (i.e., 
lim^oo p(x, t) — P q (x j) of this FP equation such that the 
inhomogeneous linear and the homogeneous nonlinear FP 
equations become unified. 



The condition 



dp{x,t) 



at 



= implies 



8F(x)p(x, 00) D d 2 

dx = 2 ft^fofM*. 00 )]- ( 5 ) 

By assuming appropriate boundary conditions (basically 
p(±oo,oo) = 0), changing variables and developing, we 
obtain 

d [g{V) e p{V) l +*] _ 2p(V) 



dV D ' 

where p{V) = p(V(x), 00). It follows that 

dp(V) 2g(Vy e p(V) 1 -i + 9BDg{V)- 1 p(V) 



dV 



D(l + V ) 



(6) 



(7) 



Now, we will first consider the 9^1 case. Let us propose 
a solution of Eq.(7) satisfying 



g{vy- e p{V)-« = C 



(8) 



with dC/dV = 0. By direct substitution in Eq. (7), we 
can easily show that p(V) = P q (V), as given by Eq. (1), 
with 



9=1 



V 



and 



/3 



2C + 9BD 
AD{l + n) 



(9) 
(10) 



p-2 



Long memory constitutes a unified mesoscopic mechanism for nonextensivity 



where 



C = [AA^ -1 ] 1 



(11) 



These results reproduce, for 9 = 0, the homogeneous 
nonlinear case discussed in [HEJE], i.e., q — 1 — rj and 
Nq~ q (3 = 2/D(l + rj), which leads to normal diffusion for 
r\ = 0. 

Eq. (9) shows that, for rj ^ 0, the limit 9 — > 1 corre- 
sponds to a singularity in the \q\ value. However, in this 
limit, we still find g-cxponcntials as stationary solutions 
of the above FP equation, in two different situations: 

a) For r] considered as a function of 9 with the leading 
term given by rj ~ a(9 — l) s , the analytical extension of 
the solution given by Eqs. (8-11) results in g-exponentials 
presenting the indices q ~ 1 for S > 1 and q = 1 + a for 
(5 = 1, with P = (2 + BD) /AD in both cases. 

b) An isolated solution (which is not a particular case 
of the previous ones) can also be found by setting rj = 
and 9 = 1 in Eq. (7). By so doing we obtain 



dp(V) 
dV 



BD]g(V)- 1 p(V) 



D 



Now, imposing the relation 



9 {v) P {vy- 1 = c. 



(12) 



(13) 



where C is independent of the potential V. Following 
the same procedure used to obtain Eqs. (9) and (10), we 
finally have p(V) = P q (V), with 



and 



(1 + BD) 
(1 + BD/2) ' 



(2 + BD) 
AD 



(14) 



(15) 



It should be noticed that this case {{r/,9) — (0, 1)) re- 
covers previous results already obtained in [TJ|3]. 
As final remarks, we note that: 

(i) For A 0, we may take A = 1 without loss of 
generality, if the (B, D) parameters are properly rescaled; 

(ii) It is known [13], that the g-Gaussians, which emerge 
for V(x) tx x 2 are normalized only for q < 3, which implies 
the following restrictions to the values of the parameters of 
the FP equation: (a) For 6 ^ 1, ^ < 2, (b) For 6 1, 
a < 2 if S = 1 and BD > -2 or BD < -4. For each 
V(x), an analysis of integrability of P q must be performed, 
in order to establish the accepted range of values of the 
parameters {rj, 9, B, D, q) (e.g., if V(x) oc \x\ p , in the limit 
|af| — ► oo, then it must be q < 1 + p); 

(iii) We stress that the Eq. (3) used in this paper cor- 
responds to the ltd form of a generalized FP equation (a 
Stratonovich approach has already been examined |10) . 
and the existence of ^-exponential distributions as station- 
ary solutions has also been proved). 

Let us emphasize that Eq. (9), here presented for the 
first time, enables us to analyze in an unique way the 



q= 1 + T|/(e-1 ) 

il=0.5 

il=1.0 

il=2.0 



Fig. 1: 



= 1 + ^3- 

1 ~ 9-1 



of rj (the upper bound q 

„2 



function of 9 for selected values 
3 herein shown corresponds to a 
potential V(x) oc x 2 , i.e., p — 2). For the special point (6* —¥ 1, 
r; — > 0), there are three different solutions for q, namely 5 = 1, 
q = 1 + a, and q = (1 + BD)/(\ + BD/2) (see the text). 

present general inhomogeneous nonlinear case, which con- 
tains, as particular cases, several FP equations used to de- 
scribe a large class of nongaussian natural systems. This 
unification constitutes the main result of this paper. 

To summarize, the nonlinear inhomogeneous FP pro- 
cess given by Eq. (3) has, as a stationary solution for any 
confining potential V(x), a probability distribution given 
by the ^-exponential P q (V), with finite values of q and /3 
(see Fig. 1 and Table 1). These results exhibit that we 
may retain long memory (basically, the probability dis- 
tribution longstandingly maintains memory of its form at 
t = 0) as a unified mesoscopic mechanism consistent with 
nonextensive statistical mechanics. 

We acknowledge useful remarks by A.O. Sousa, as well 
as partial financial support by Fapern, Faperj and CNPq 
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FOKKER 
-PLANCK 



Linear 

(ri->0) 



Nonlinear 

(ry^O) 



Homogeneous 

(9 = 0) 




q = l- T] 



Inhomogeneous 

(MO) 



If r] ~ a(6 - l) s and 6 - 
q = 1 for S > 1 
q = 1 + a for 5 = 1 



1- 



9 = 1 + 



If j] is strictly and 

_ 2(1 + BD) 
— 2+BD 



Table 1: Nonlinear inhomogeneous Fokker-Planck equation with q-exponential stationary-state distribution. The present 9 = 
result recovers that of Eqs. (2.8) and (2.9) of [4] with the notation correspondence (F, Q, D, K) in 4] (P, D, N q , F) here. It 
also recovers that of [5] by using there /i = 1 and v = 1 + r\. The result presented by the isolated solution at (8, rj) = (1,0) 
recovers that cited (for the Ito approach) in [T] with the notation correspondence (— ) in [lj (^p) here. A comparison is 
also possible through Eq. (11) of [3] by using the notation correspondence (K, D,U(x)) in [3] -f-J- (F, Dg(V) , V(x)) here. 
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